Tunneling lifetime of trapped condensates 
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Abstract 



Using the complex scaling transformation we calculate the tunneling lifetime 
of a condensate inside a magnetic/optical trap. We show that by varying 

o. 

the scattering length, the external potential acts like a 'selective membrane' 



which controls the direction of the flux of the cold atoms through the barriers 



and thereby controls the size of the stable condensate inside the trap. 
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Tunneling of Bose-Einstein condensates (BEC) through a dynamical potential barrier has 
been demonstrated in the experiments of Phillips' group [1]. Recently it has been shown that 
the tunneling of cold atoms can be controlled from full suppression to strong enhancement 
by varying slightly the experimental parameters [2]. In contrast to the previous studies 
where the tunneling is through a dynamical potential barrier, we discuss here the tunneling 
of the BEC from a potential well to the continuum through a static potential barrier. In all 
BEC experiments the condensates are trapped inside an open external potential well which 
is embedded in between potential barriers. Therefore, in the case of a positive scattering 
length (repulsive atom-atom interactions) a given external open potential well cannot trap a 
condensate with more than a finite critical number of atoms. As the number of the atoms in 
the condensate, N, is increased, fractions of the condensate can tunnel through the potential 
barriers. 

What is the lifetime of a trapped condensate as N is increased ? Can we control the 
number of atoms in the condensate by varying slightly the scattering length ? Using the 
complex scaling transformation which is based on rigorous mathematical ground [3,4] and 
has been used before for calculating the lifetimes of metastable states in atomic, molecular 
and nuclear physics [5,6] we show that the lifetime of the condensates can be calculated and 
the answer to the other question is positive as well. 

For the sake of clarity and without loss of generality let us represent a simple one- 
dimensional external potential, V ext (x) = (x 2 /2 — 0.8) exp(— O.lrr 2 ). For a single atom this 
external potential supports only a single bound state. Consequently, in the absence of 
interaction between the atoms of the BEC, the condensate can be found permanently trapped 
inside the potential well or be in one of the metastable states (so called resonances) where 
the atoms can tunnel out of the trap (see inset in Fig. 1). We shall see that in the presence 
of repulsive interaction between the atoms the situation changes. Whether the condensate 
can be in a bound state at all depends on the interaction strength and on the number of 
atoms. 

Before discussing our results we discuss in the following metastable states and how they 



are determined by complex scaling in the case of condensates. In hermitian quantum mechan- 
ics metastable states are associated with propagation of wavepackets escaping the potential 
well. However, in quantum scattering theory it is possible to associate metastable states 
with the solutions of the time-independent Schrodinger equation. The metastable states, as 
are the bound states, are associated with outgoing boundary conditions. In the case of BEC 
the metastable states are embedded in the continuum part of the external potential. It is 
known from scattering theory that outgoing solutions which are embedded in the contin- 
uum are obtained only when the wave vector, k, gets complex values k = \k\ exp(— ia), and 
hence the eigenvalues of the Hamiltonian, {£}, also become complex quantities [9]. Complex 
eigenvalues are obtained since the corresponding eigenfunctions diverge exponentially, i.e., 
^res ~^ exp(±i|/c| cos(a)a;) exp(±|/c| sin(a)rr) — ► +00 as x — »■ ±00, and are not embedded 
in the hermitian domain of the Hamiltonian. {ReS} are associated with the energies of 
the metastable system and {—2ImS}, are associated with the rates of the decay (inverse 
lifetimes). By complex scaling of the internal coordinates in the Hamiltonian (in our case 
x — > xexp(i9) where 9 > a) the metastable states are taken back into the Hilbert space 
and become square integrable (like the bound states). This approach enabled developing 
the quantum theory for non- hermitian Hamiltonians (see the Reviews in Ref. [5]) and Ref. 
[6]. 

How to compute the resonances of BEC ? The BEC consists of N atoms that are 
assumed to have contact potential interactions, Uq8{tj — ry), where Uq = Anctoh 2 jM 
and ao is the s-wave scattering length [7]. To calculate the bound and the resonance 
(metastable) states, as explained above, we carry out the following similarity transformation 
[6], Hbec = SHbecS 1 , where S is the complex scaling operator. The result of the com- 
plex scaling transformation is, Hbec(&) — —^ 2 /2MJ2j ( ex P( — 2i9)V 2 + V ext (exp(i9)fj)) + 



Uoexp(—i8) J2j,j'>j ${rj — fji). The bound and the resonance states are obtained by solving 
the complex eigenvalue problem, H.bec(Q)^(9) — £^(#)> where ^(9) are square-integrable 
functions and £ are 6- independent eigenvalues. The bound states (i.e., ImS = 0) are ob- 
tained in our case for any value of < 9 < ir/4. The resonance energies and widths per 



atom, S/N = E — i/2T, are obtained when 7r/4 > 9 > 9 C . The value of the critical scaling 
angle is given by 29 c = arctan2i?/T [5,6]. Let us return to the statement that ^(9) are 
square-integrable functions. As we will see later this is an important point in the derivation 
of the non-linear Schrodinger equation for complex scaled Hamiltonians. Square integrabil- 
ity implies that the corresponding functions decay exponentially to zero. But this is not 
enough. We need to define the inner product. In hermitian quantum mechanics where 
W = H we use the common scalar product. In the case of complex scaled Hamiltonians 
W ^ H. In our case (as in most studied physical cases for complex scaled Hamiltonians [6] 
), Ti-BEci®) = Wbec(Q)i an d therefore the c-product rather than the scalar product should 
be used [6] such that, (^(0)1^(0)) = (**, (0)|*,,(0)) = S fJ , where V r (0) and ^(0) are 
any bound or resonance eigenfunctions of 1~Lbec{Q) as defined above. Note that here we 
assume that for = (i.e., hermitian quantum mechanics) the eigenfunctions of H are real. 
The non-linear Schrodinger equation, also known as the Gross-Pitiavesky equation, is 
commonly applied to condensates [7]. It is a mean field approximation to the full many-body 
Hamiltonian of the problem as the Hartree-Fock approximation is for fermions. All atoms 
are assumed to occupy the same orbital and the total wavefunction is a product of these 
orbitals [7]. Therefore, S =< V*(6)\Hbec\V(0) >, where V(0) = I\f=i M^j) ■ Minimizing 
£ with respect to the orbital, keeping in mind the definition of the inner product, leads to 
the complex scaled version of the non-linear Schrodinger equation: 

(H (9) + Uexp(-iO)tf) h = {^- \l)<\>e, (1) 

where, H (9) = -h 2 exp(-2i0)/2MV 2 + V ext (exp(i9)r) and U = U (N-1) = (Aixa Q h 2 )(N - 
1)/2M is the relevant non-linear parameter. The complex chemical potential of the 
metastable state, \x — |7 (7 = for bound states), is associated with the complex energy of 
the BEC per atom, S/N = E — i/2T, as given by, 

E- t -T = f,- t - 1 + ^ exp(-t9) J <f)tdr. (2) 

Z Z Z J all— space 

Within the framework of the Gross-Pitiavesky approximation resonances can only be 



calculated when the complex scaling transformation is applied to the BEC N-body Hamil- 
tonian before applying the mean field approximation. To see why one cannot apply first 
the mean-field approximation and only later the complex scaling transformation, we start 
from the usual Gross-Pitiavesky Hamiltonian and scale it. The result is given by, Ti.ap(0) = 
SHgpS- 1 = -exp(-2i6)hy2MV 2 + V ext (exp(-i6)r) + U[S<j)* r JlS<j) res }. The resonance solu- 
tion of Eq. 1 is square integrable. It is, however, associated with the exponentially diverging 
unsealed resonance solution, <p res , such that, (pe = S(p res — ► exp(+i|/c| exp(+i(9 — a))r) where 
9 > a. Therefore, [S(f>* es ] [S<f> res ] = exp(— i2\k\ sin(0) sin(ai)r) exp(+2|/c| cos(#) sin(ai)r) — > 
+00 as r — > 00. Namely, the complex scaled non-linear potential term U[S<j)* es }[S<j) res ] ap- 
pearing in TicpiQ) diverges as r — > 00 (although <p 9 decays exponentially to zero) and its 
meaning is unclear. For U > all states become bound and for U < the potential is 
unbound from below. Therefore, it is impossible by this approach to get an eigenfunction 
which is associated with the complex scaled resonance eigenfunction. Note that this diffi- 
culty has been avoided in a previous study of bound-resonance state transitions using the 
Gross-Pitaevski equation by inserting complex absorbing potentials in the regions where it 
has been assumed that the atoms do not interact with one another [8]. However, here we 
do not use this approximation and the atoms can interact everywhere and not only inside 
the external potential well. 

We solved Eq. 1 for the ID potential defined above, by increasing adiabatically the 
non-linear parameter U. We used 400 particle-in-a-box basis functions as a basis set with 
the box size of L = 50a.u.. In each step of the calculations we carried out self-consistent 
field iterative calculations to get converged results in 8 significant figures where the scaling 
angle 9 has been varied from 0.3rad to 0.7rad. By solving Eq. 1 we calculated the complex 
chemical-potential, /i — i/27 and by solving Eq. 2 we computated the complex-energy per 
atom E/N = E- i/2T. 

In Fig. 1 we show T and 7 as functions of U. We mark by arrow the critical value 
U c = 0.8279 for which the bound-resonance transition occurs. For U < U c the system is 
bound, i.e. trapped forever in the well, and for U > U c it is metastable and tunneling 



through the barriers occurs. The trap potential is depicted in the inset of the figure. The 
corresponding real parts fi(U) of the chemical potential and E(U) of the energy are shown 
in Fig. 2 as functions of U. As one can see, at the critical value Uc where the resonances are 
"born" the chemical potential vanishes, /i c = fi(U c ) = 0. This is an expected result. We may 
interpret fi as the energy needed to take a single atom of the condensate out of the potential 
well. 7 is the corresponding rate of decay, i.e. I/7 is the tunneling time of a single atom 
with a chemical potential \x. The lifetime of the BEC resonance state is l/(iVT), which can 
be short if the condensate possesses many atoms. In a condensate all atoms are equivalent 
and each of them can tunnel with the same probability. Interestingly, T > 7 for all values of 
the non-linear coupling U > U c as can be seen in Fig. 1. The reason is persumably that 1/T 
is the lifetime per particle, and each of the BEC's particles can tunnel through the barrier 
in different ways, as a single particle or together with several other equivalent particles. As 
mentioned above, I/7, on the other hand, is associated with the tunneling of a single atom 
only. 

How to visualize the decay of a many-body system like a BEC? Once the BEC is in 
a metastable state, a fraction of the atoms consisting of (l-X)N atoms tunnels out of the 
potential well into the continuum. This leads to a stabilization of the system where XN 
atoms remain in the trap. Since all the atoms in the trap repel each other, the energy of 
the XN atoms is lower than that of the N atoms. Tunneling proceeds until the fraction of 
atoms which remains in the trap forms a bound state. The number of atoms which will 
tunnel clearly depends on the energy of the condensate. However, within the framework of 
the Gross-Pitaevski mean- field approximation the non-linear parameter is U — Uq(N — 1) 
where N is the total number of atoms, and one cannot tell at a given value of U which 
fraction of N will tunnel. As we will show below we can enlighten this problem by taking 
into consideration the idea that not always the best mean-field for condensates is obtained 
for the case where all identical bosons of a condensate reside in a single orbital [10]. 

We consider here the scenario where the ground state \1/ is a product of two types of spatial 
orbitals. There are ri\ atoms which occupy the orbital and 112 = N — n\ atoms occupy the 
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X orbital. We associate the n\ atoms in the orbital <fi with the fraction of the condensate 
X = rii/N which are located in the well and the n 2 atoms in the orbital x with the fraction 
(1 — X) = n-ij 'N which has tunneled through the barriers into the continuum. Consequently, 
is located in the trap and x we U outside the trap where V ex t = 0, and we may consider 
these two orbitals to not overlap. Since U = Uq(N — 1) ~ UqN, the energy per atom in the 
trap E^,(X,U) is nothing but E(XU) which is reported in Fig. 2. Analogously, the energy 
per atom in the continuum is given by E x ((l — X) , U). The resulting energy per atom of the 
whole system at a given U and X is thus E BEC (X, U) = XE(XU) + (1-X)E X ((1-X), U). 
Because our condensate is repulsive, the energy E x is known to be equal to zero [7] and we 
obtain the final result, 

E BEC (X,U)=XE(XU), (3) 

which makes clear that Ebec{X, U) can be derived from the curve E(U) in Fig 2. In complete 
analogy the rate of decay into the continuum per atom corresponding to E BEC (X, U) is given 
by XT(XU) and is also already computed, see Fig. 1. 

While for a given value of U the Gross-Pitiavesky energy E(U) is just a number and 
does not provide us with the knowledge on how many atoms have tunneled, Ebec{X, U) is 
the key to this information. In Fig. 3 we show the energy per atom of the condensate as a 
function of X for different values of U. The curves at the bottom of the figure are for small 
values of U, those at the top for larger values of U. Each of the curves exhibits a minimum at 
X C (U) and these minima play a central role in the understanding of the tunneling process. 
These minima are marked by solid dots. Let us consider a single curve in Fig. 3 for which 
X c is smaller than 1 (the value X=l is marked in the figure by a vertical line). X=l implies 
that we have put N atoms in the condensate and is thus our starting point. Because of the 
variational principle, the condensate will minimize its energy by letting a fraction 1 — X c 
of its atoms tunnel into the continuum and keep the fraction X c in the well at which Ebec 
takes on its minimum. 

How does this appealing picture relate to the decay rate XT(XU) ? We know already 
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from the scaled mean-field approach that T(U) changes from being zero to non-zero at the 
bound-resonance transition point U c . Therefore, for a given value of U we find a particular 
value of X such that XU = U c , and this X tells us at which fraction of the condensate the 
system is just still bound. If this is the case, there should be an intimate relation between 
this particular value of X and X c at which Ebec has its minimum. Indeed, we find that both 
values of X are identical, i.e. X c = U c /U. Using the value U c = 0.8279 found above, one 
readily reproduces the values of X at which any curve Ebec(X,U) takes on its minimum 
for a given U. In particular, Ebec{X,U c ) -which is the blue curve in Fig. 3- exhibits its 
minimum at X c = 1. This value implies that for U c the condensate with N atoms is bound, 
while for U > U c the fraction 1 — X c tunnels to make the remaining fraction X c < 1 bound. 

It is essential to note that the critical values X c and U c can be determined from the 
above analysis without using the complex scaled mean- field results. The curves Ebec{X, U) 
shown in Fig. 3 can be computed via Eq.3 for all values of U from X=0 up to X c using bound 
state calculations only. This is a success of the two-orbital picture [10] used above to derive 
Eq. 3. 

Once U > U c , increasing the s-wave scattering length ao, for instance, by applying an 
external magnetic field to adjust the relative energy of different internal states of the atoms 
[11], leads to an increase of U and hence to a reduction of the number of atoms inside 
the potential well (see Fig. 3). On the other hand, if we decrease the scattering length 
ao, the trap can accommodate more atoms. Consequently, if there is a reservoir of cold 
atoms outside the trap, some of them can tunnel through the barriers into the trap thus 
increasing the number of atoms inside the trap. Note that in Fig. 3 the minima of Ebec f° r 
U < U c are at X c > 1, i.e, the condensate inside the trap is further stabilized if atoms are 
added. In this way the size of the trapped condensate can be controlled and the trap acts as a 
"controllable membrane" by varying the s-wave scattering length (or by varying the depth of 
the trap potential). We hope that these fascinating results will stimulate new experiments. 
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FIG. 1. The rate of decay 7 of a single atom and the total rate of decay per atom T as a 
function of the non-linear parameter U (see Eq. 1 and text). The inset shows the external trap 
potential used. 
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FIG. 2. The chemical potential //(the real part of the complex eigenvalue in Eq. 1) and the 
mean-field energy of the BEC per atom E (the real part of the complex energy S/N, see Eq. 2 in 
the text) as a function of the non-linear parameter U. 
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FIG. 3. The energy per atom Ebec(X,U) as a function of the fraction of atoms X which 
remains in the trap while a fraction 1-X of the condensate has tunneled through the barriers into 
the continuum. Each curve shown is for a different value of the non-linear parameter U. From 
bottom to top the j-th curve in black is associated with U=0.5 + 0.05(j-l). The blue curve is for 
U = U c . The minima of the Ebec curves are at X=Xc (solid dots) and play a central role in 
understanding the tunneling (see text). 
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